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NbP is one member of a new class of nodal loop semimetals characterized by the cooperative 
effects of spin-orbit coupling (SOC) and a lack of inversion center. Here transport and spectroscopic 
properties of NbP are evaluated using density functional theory methods. SOC together with the 
lack of inversion symmetry splits degeneracies, giving rise to “Russian doll nested” Fermi surfaces 
containing 4x10“'* electron (hole) carriers/f.u. Due to the modest SOC strength in Nb, the Fermi 
surfaces map out the Weyl nodal loops. Calculated structure around T*«100 K in transport prop¬ 
erties reproduces well the observed transport behavior only when SOC is included, attesting to the 
precision of the (delicate) calculations and the stoichiometry of the samples. Low energy collective 
electron-hole excitations (plasmons) in the 20-60 meV range result from the nodal loop splitting. 

PACS numbers: 71.20.Be,71.18.-l-y,78.20.Ci,72.15.Eb 


I. INTRODUCTION 


Weyl semimetals (WSs) have Weyl points, comprised 
of a linear crossing of valence and conduction bands [l| 
near the Fermi energy E^r. WSs can be realized from 
a lack of time-reversal symmetry (magnetic ordering, or 
applied field B), GHi or from no inversion symmetry in 
three-dimensional compounds, a a and show peculiar 
features, such as Fermi arcs from surface states connect¬ 
ing two Weyl points with different chirality. aai Ex¬ 
perimental evidence has been obtained using x-ray angle- 



FIG. 1: (Color online) Body-centered tetragonal liimd struc¬ 
ture of NbP. It consists of edge-sharing NbPe trigonal prisms, 
which are 90° rotated in the next layer. The lack of inversion 
symmetry is evident: the prisms all “point” upward. 
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resolved photoemission spectroscopy. @ Extremely large 
magnetoresistivity (XMR) also occurs in some WSs, with 
a recent analysis by Burkov Q concluding that XMR 
should be expected in high mobility semimetals where 
Weyl or Dirac points lie near the Fermi surface (FS). 
XMR was observed also in the semimetal WTe 2 , show¬ 
ing unsaturated MR up to 60 T,0 attributed to precise 
electron and hole carrier compensation, as in Bi and other 
stoichiometric semimetals. 

Recently a “nodal loop” Weyl semimetal phase with 
potentially distinctive properties has been proposed 
on theoretical grounds [J, in the nonmagnetic non¬ 
centrosymmetric isovalent class TaAs, TaP, NbAs, and 
NbP, whose structure is shown in Fig. 1. A previ¬ 
ous proposal of nodal loop WSs by Phillms and Aji in¬ 
voked time-reversal symmetry breaking[8j rather than 
the lack of inversion symmetry of this class of pnictides. 
Since these suggestions, XMR has been observed in both 
TaAs[^ [Toll and NbP.flO, [HI NbP, upon which we focus 
here, shows XMR of 8,500 at 2 K for i? = 9 T, increas¬ 
ing to 3 X 10* at 4 K at B=30 T, to 10® at 60 T, and 
remains unsaturated. [HI These colossal XMR values re¬ 
flect the excellent conduction at low temperatures, with 
resistivity as low as 0.1/iD cm, and a metal to insulator 
crossover with field at 100 K and below. A high mobility 
/i=5 X 10® cm^/V-s and low carrier density of 1.5 x 10*® 
cm“® at low T was inferred from the conventional single 
band expressions, which however are not quantitative for 
multiband semimetals. As temperature increases, both 
the Hall coefficient Rh changes sign from electron-like to 
hole-like and the mobility rapidly changes in the 50-150 K 
range, correlating with the metal to insulator crossover. 
A multiband model and additional data will be neces¬ 
sary to sort out individual carrier densities and mobil¬ 
ities. NbP has the smallest spin-orbit coupling (SOC) 
of this class, making it the most delicate example of the 
features of nodal loop WSs. 

From both theoretical and experimental viewpoints. 










nodal loop semimetals comprise a new class of materials 
with topological behavior. In this paper, we focus on ef¬ 
fects of SOC on their electronic properties and evaluate 
the optical and transport coefficients of NbP, using den¬ 
sity functional theory based methods that enable a direct 
theory-experiment comparison. A central feature that 
emerges is electronic fine structure driven by the combi¬ 
nation of SOC and lack of an inversion center, a feature 
elaborated on by Samokhin[l^ from a model viewpoint. 
This splitting in a seinimetal with small closed Fermi sur¬ 
faces (FSs) results in a “Russian doll nested” pair of FSs, 
and new low energy interband transition strength pilfered 
from the Drude strength by non-centrosymmetry. The 
FS character is clarified and low energy collective elec¬ 
tronic excitations (plasmons) arising from the electronic 
fine structure are predicted. While plasmons have been 
studied in conventional WSs. [l3l - [l^ their appearance in 
nodal loop WSs as well as several other properties dis¬ 
cussed here are distinctive when the special “nodal loop 
semimetal” character is encountered. 


II. STRUCTURE AND THEORETICAL 
METHODS 


As shown in Fig. [I] NbP has a body-centered tetrag¬ 
onal lattice of the space group lAimd (^109). [a This 
space group has 4mm point group, is non-symmorphic, 
and (importantly) lacks inversion symmetry. The struc¬ 
ture can be pictured as columns of face-shared NbPg 
trigonal prisms oriented along the a-axis, edge-shared 
along the &-axis in one layer. The trigonal columns 
are 90° rotated in the next layer along the c-axis, giv¬ 
ing overall tetragonal symmetry. In the structure with 
the experiment lattice parameters [l7i of a = 3.3324 A, 
c = 11.3705 A, both Nb and P ions he at 4a posi¬ 
tions (0,0, z): ZNb = 0, zp=0.4176; the site symmetry 
is 2mm. Our optimized positions using the generalized 
gradient approximation (GGA) with the Perdew-Burke- 
Ernzerhof (PBE) exchange-correlation functional (l8j| in 
FPLo[l^ coincided with the experimental values. This 
structure leads to nearly identical Nb-P bond lengths of 
2.53±0.01 A and P-Nb-P bond angles around 82°. 

Calculations with the experimental structure were 
performed using the all-electron full-potential code 
wien2k,[2^ with selected results confirmed with fplo. 
All results are based on the PBE-GGA exchange- 
correlation functional. One objective is to determine the 
combined effects of SOC and the lack of inversion center, 
since delicate features around Ep are sensitive to SOC 
and become important. 

Calculation of optical properties including SOC is 
available in wien2k.[21| The dielectric function eij(uj), 
with only diagonal Coq = thb and tec components due to 
tetragonal symmetry, is decomposed into the intra- and 
inter-band contributions. The imaginary parts of each 



FIG. 2: (Color online) (a) GGA bands in the range of -1.5 
eV to 1.5 eV and (b) GGA-I-SOC bands in the range of -0.2 
eV to 0.2 eV, with Ep set to zero. The symmetry points 
can be identified in Fig. [S] The upper right panel enlarge¬ 
ment shows the splitting by ~10 meV of bands by SOC in 
this non-centrosymmetric structure, (c) Density of states 
(GGA-I-SOC) showing the pseudogap centered 20 meV be¬ 
low Ep. 


contribution are given byilil 
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where p is the momentum operator, and £c^ and are 
energies of the occupied Vj: and unoccupied Cj^ orbitals, 
respectively. In the intraband term, which contains the 
Drude divergence for w —>■ 0, ^pjj is the Drude plasma 
frequency and F (chosen to be 10 meV here) is an inverse 
scattering lifetime r. The corresponding real parts are 
obtained by the Kramers-Kronig relation. 

The transport calculations were carried out using 
semiclassical Bloch-Boltzmann transport theory, imple¬ 
mented in the BoltzTraP code [13 that is interfaced to 
WIEN2k. No approximations beyond the constant scat¬ 
tering time approximation are made. In all calculations 
performed here, the Brillouin zone was sampled with a 
dense fc-mesh up to 60 x 60 x 60 to treat the semimetallic 
fine structure carefully. In wien2k, the basis size was de¬ 
termined by RmtKmax = 8, and augmented plane-wave 
sphere radii of 2.5 Bohr (Nb) and 2.14 Bohr (P) were 
used. 


HI. RESULTS 


A. Electronic fine structure 

We first address the electronic structure obtained with 
fine fc-point and energy meshes. Figure [2] shows the GGA 
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FIG. 3: (Color online) NbP Fermi surfaces (FSs) including 
SOC and shown from two angles, the lower one viewed along 
the kx axis to reveal the thinness of the surfaces in the mirror 
planes. The extended zone representation makes the toroidal 
origin of the FSs evident. The hole (blue) boomerang FSs 
form bent and ‘sharply’ tipped boomerang surfaces, the elec¬ 
tron (yellow) new moon surfaces are similar but extend over 
a larger angle of the torus. 


and GGA-I-SOC band structures near £^7-. Semimetallic 
character with a band crossing along at least three sym¬ 
metry directions appears, with valence-conduction over¬ 
laps mostly of the order of 150 meV. This band behavior 
leads to a pseudogap centered 20 meV below Ep (Fig. [2|), 
with a very small density of states (DOS) N{Ef)=0-04:6 
states/eV per f.u. While these results are similar to the 
previous report, ami we examine the origins and conse¬ 
quences more completely. 

Inclusion of SOG leads to unconventional conse¬ 
quences. As usual, SOC converts many band crossings 
along symmetry lines into anticrossings; here gapping 
bands by a few tens of meV near EiT-. However, the lack 


of inversion symmetry leads to lifting of band degener¬ 
acy over entire symmetry planes, a with splittings in the 
20-200 meV range. In spite of the modest magnitude of 
SOC in 4(i Nb, SOC will play an important role in deter¬ 
mining both the thermal and spectroscopic properties of 
NbP. 


B. Fermi surfaces 

The FSs obtained from GGA-bSOC are displayed in 
Fig. [31 see the caption for a description. They are sim¬ 
ilar to those presented earlier, but they are displayed 
in an extended zone manner to illustrate their toroidal 
origin. Weng et al. have described in detail^ the un¬ 
derlying features of the Weyl nature of NbP and isova- 
lent partner compounds. Without SOC, mirror symme¬ 
try protected nodal loops lie precisely in the mirror planes 
kx = 0 and ky=0. The nodal loops consist of band de¬ 
generacies Ei^sik) = e 2 ,s{k) (s is spin degeneracy). The 
locus of k points at this intersection forms the ring in 
the PZN mirror plane, not at constant energy but, for 
these compounds, in a small energy range including Ep- 
The FS intersects the nodal loop, so it must include an 
even number of doubly degenerate Fermi points where 
Si,s{k) = £ 2 ,s{k) = Ep. Such nodal loops have also been 
found in two dimensional SrVOa nanolayers. [2^ 

In this non-centrosymmetric structure, SOC splits the 
bands forming the nodal loops, leaving two symmetry- 
inequivalent pairs of Weyl points (12 in all) near but 
not on the FSs, the points being provided explicitly by 
Weng et al.Q One lies in the kz = 0 plane, the other 
at a general k point lying near the PNZ plane. The 
band splitting is highly anisotropic; each FS (Fig. [3]) 
consists of a flattened torus (annulus) pinched off into 
one “boomerang” containing holes and one “new moon” 
containing electrons. Small doping will alter the position 
of the pinching, hence changing the number of electron 
and hole carriers. 

For a stoichiometric compound the FSs contain an 
equal number of holes and electrons, and from calcu¬ 
lated band fillings we obtain 4.2x10“^ electron (hole) 
carriers per f.u. This value corresponds to high veloc¬ 
ity carriers of each sign separated on average by 13 lat¬ 
tice constants in each direction. The rms Fermi veloc¬ 
ities are ?^F.aa=3.7, r'F,cc=l-6, in 10^ cm/s, a factor of 
2.35 in anisotropy. Note that the velocities are typical of 
metallic compounds, not semimetals, due to the disper¬ 
sive bands crossing Ef, i-e. the Weyl character (anal¬ 
ogous to graphene). For the Drude plasma energy, the 
semimetallic value of N{Ep) but normal metal velocities 
lead to Dp_QQ=1.0 eV, Dp^cc=0.4 eV. 

What cannot be seen in Fig. |3| is the doubling of 
the FSs resulting from the interplay of SOC with the 
non-centrosymmetric crystalline symmetry. With the 
band degeneracy split, Fermi surfaces become “Rus¬ 
sian doll nested” pairs, one inside the other. This 


3 













T[K] 

FIG. 4: (Color online) (a) Temperature (T) dependent Hall 
coefficients Rh{T) calculated in both GGA and GGA+SOC, 
which show similar T-dependence but distinctions in the zero 
crossing temperature; T=125 K in GGA+SOC, compared to 
T=230 K in GGA. (b) T-dependent total density of states 
N{T,iJ,) at chemical potential N{T,jj,) shows a min¬ 

imum around T=100 K in GGA+SOC, while that of GGA 
increases as increasing T. (c) with respect to Ep at T=0 
K. 


band-splitting effect has been discussed for Pt-based 
superconductors [ 2 ^ and for electronic properties more 
generally, [l^ The band spitting aspect is effective both 
in fc-space (closely nested Fermi surfaces) and in energy, 
with consequences discussed below. 


C. Transport properties 

Transport behavior in semimetals is very sensitive to 
Fermiology and to regions of high curvature, but we 
have found that calculated values based on quasiclassical 
Bloch-Boltzmaim theory agree very well with observed 
behavior. m The Hall tensor components Rh(T) includ¬ 
ing SOC, displayed in Fig. lUJa), show a change in sign at 
rc=125 K, precisely where the experimental R// changes 
sign. [iTj l This change of sign is a strong reminder that in a 
multiband, compensated semimetal, R// bears no direct 
relation to carrier densities. This statement is partic¬ 
ularly true when the FSs are strongly non-ellipsoidal and 
extremely anisotropic. 

The experiment-theory agreement reflects not only the 
precision of the calculations but also the excellent stoi¬ 
chiometry of the samples; the stoichiometry could not be 
established precisely but is supported also by the high 
mobility. The magnitude of Rh{T) around 250 K is 
within a factor of 2 of the experimental value. This is 




FIG. 5: (Golor online) Temperature and polarization depen¬ 
dent (a) conductivity tensor over scattering time, in 10^'^ (fl- 
m-s)“^, and (b) Seebeck coefficients Su [/rV/cm], note the 
minimum in the latter in the 100-150 K range, (c) The See¬ 
beck coefficient 5'aa(p; T) versus band filling /r at various tem¬ 
peratures, within GGA+SOG. The change at very small levels 
of doping is extremely rapid, reflecting the movement away 
from charge compensation; Sec behaves in the same manner. 
The temperature variation is monotonic. 


really excellent agreement, because of the compensating 
contributions from electron and hole surfaces, which have 
regions of large curvature. The expression for Rh in 
Bloch-Boltzmann theory has the interpretation of giving 
it as the average of the curvature over the FSs. A strong 
temperature dependence and a likely change of sign with 
temperature should be anticipated in nodal loop WSs. 
The experimental Ru becomes at the lowest tempera¬ 
ture two orders of magnitude larger in magnitude than 
our calculated value, and it is clear from the “pointed” 
FSs in Fig. [3] that evaluating Rh (which as mentioned is 
an average of the FS curvature) at T=0 is a substantial 
numerical challenge. 

The chemical potential /i(T), displayed in Fig. SKc), 
decreases with increasing T reflecting the strong differ¬ 
ences in the electron and hole densities of states. This 
variation is different in sign from another XMR com¬ 
pound WTe9.[^ reflecting different DOSs of the valence 
and conduction bands of the two compounds. The DOS 
at the chemical potential shown in Fig. HKb), is re¬ 

duced and achieves a minimum around 100 K before turn¬ 
ing around. Figure [SKa) displays the conductivity tensor 
components CT^i/r, which incorporates the anisotropy of 
Dpjj, which is ~5 times larger in-plane. 

The T-dependent Seebeck coefficients, shown in Fig. 
[SKb), have net n-type sign. They peak (in magnitude) 
at 100-150 K before again becoming small. The depen- 
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FIG. 6: (Color online) The polarization-dependent energy loss 
function for NbP assuming a scattering rate F = ?i/r of 10 
meV. Inset: the c-axis plasmon at 20 meV bifurcates into a 
pair of plasmons at 40 meV and 60 meV for in-plane polariza¬ 
tion. Above 0.2 eV the loss function is not greatly anisotropic 
and is insensitive to F. 


dence on chemical potential which can be varied by 
doping or by gating, is extremely strong near stoichiom¬ 
etry [see Fig. 4(c)]. The near-perfect cancellation of 
electron and hole contributions (it is perfect Sjj = 0 for 
/iss-lS meV) at stoichiometry is destroyed, with a maxi¬ 
mum of |>S'(/r)| being attained around 0.5% hole doping, 
or half of that for electron doping. The sensitivity makes 
the thermopower an important gauge of the degree of 
stoichiometry of samples, which otherwise can be very 
difficult to determine at this low level. 


D. Optical properties 

The band splitting induced by SOC and the lack of in¬ 
version symmetry have qualitative consequences for the 
optical properties. In conventionalWSs where Fermi sur¬ 
faces are simple and model band structures work well, the 
dielectric response and plasmon signatures of chirality 
have been well studied, [isl-llq Nodal loop WSs bring new 
processes to the low-energy, small-g dielectric response. 
They are compensated semimetals thus providing itiner¬ 
ant electron and hole gases. They are qualitatively differ¬ 
ent from other semimetals (including conventional WSs) 
by having highly non-parabolic and extended Fermi sur¬ 
faces, and of course by having very closely nested Fermi 
surfaces. New dielectric phenomena must emerge. 

Without either SOC or inversion symmetry breaking, 
bands are doubly degenerate and contribute accordingly 
to the intraband Drude response. In nodal loop WSs 
half of this weight is transferred to finite-g, finite-w re¬ 
sponse. The changes will be appreciable only at small 
energies, and they should involve new response at very 
small q wavevectors that connect the finely nested Fermi 
surfaces, which are on the order of I0“^7r/a. The Drude 
response will be relatively unchanged apart from the fac¬ 
tor of 2 reduction in weight, but it may be smeared into 


the new small-w, small-g structures. 

The energy loss function, shown in Fig. [5] with scat¬ 
tering rate F=I0 meV, displays low energy plasmons 
when SOC is included, at 20 meV for c-axis polariza¬ 
tion and double in-plane plasmon peaks at 40 meV and 
60 meV. The strength of these peaks, and to some degree 
their positions, are dependent on the presumed scatter¬ 
ing rate F. The anisotropy we calculate at low energy, 
based on the strong anisotropy of the Fermi surfaces, 
may be consistent with the angle-dependence of the MR 
measurements. [ll| Given the reported high mobility, the 
plasmon peaks may be much sharper at low tempera¬ 
ture than calculated here. The peaks can be expected 
also to have temperature dependence, because the vary¬ 
ing /r(T) induces changes and F will be T-dependent as 
well. Above 200 meV, we do not expect significant T- 
dependence of the spectral behavior. Far-infrared spec¬ 
troscopy may reveal unusual behavior as temperature, 
field, polarization, and carrier density are varied and 
tuned. 


IV. DISCUSSION AND SUMMARY 

Our study has illuminated the origins of the unusual 
and delicate Fermi surfaces of NbP. Without considera¬ 
tion of spin-orbit coupling, the electronic structure in¬ 
cludes nodal loops 0 lying in the mirror planes that 
lead to point Fermi surfaces as well as small electron 
and hole pockets. In this non-centrosymmetric structure, 
spin-orbit coupling removes the spin degeneracy, result¬ 
ing in Russian doll nested pairs of Fermi surfaces with a 
pinched-off annulus topology, with electron surfaces con¬ 
verting to hole surfaces across the pinch. The resulting 
electronic fine structure accounts well for the change in 
sign of the Hall coefficient at 100 K, which is also the 
temperature range in which the magnetoresistance be¬ 
gins to grow in size. These features identify 10 meV as a 
relevant energy scale. 

The unusual Fermi surface also accounts for the lack 
of universal scaling of XMR in NbP and TaP; spin-orbit 
coupling is much larger in TaP and the Fermi surfaces 
are sensitive to this. Finally, the fine structure in NbP 
results in the appearance of polarization dependent low 
energy plasmons in the 20-60 meV range. This study sub¬ 
stantially broadens the understanding of how WS behav¬ 
ior impacts the physical properties of NbP, which should 
also apply in similar form to NbAs, TaP, and TaAs. 
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